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A general Landauer-Biittiker-type current formula is derived, which can be applied to the ferro- 
magnet (F)/Normal metal (N) /superconductor (S), F/N/N, N/N/S and N/N/N systems, even 
in the presence of interactions in the central region. 
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I. INTRODUCTION 



Electronic transport in mesoscopic systems or 
nanoscopic structures has received extensive, and inten- 
sive theoretical and experimental attentionEJ. In meso- 
scopic systems the sample size is smaller than the phase 
coherent length, electrons retain their phase when trav- 
elling through the sample. In the ballistic limit, i.e., 
when the dimensions of the sample is smaller than the 
mean free path, electrons can traverse the system with- 
out any scattering. In contrast to macroscopic systems, 
the conductance of mesoscopic systems is sample-specific, 
since electron wavefunctions are strongly dependent on 
the form of the boundary of the sample and the con- 
figuration of scatterers located within the sample. To 
calculate the conductance of mesoscopic samples, one 
should at first consider the wa^je nature of electrons. The 
Boltzmann transport equationo is obviously inappropri- 
ate, since the preassumption that electrons can be viewed 
as classical particles does not hold at a mesoscopic scale 
due to the Heinsenberg uncertainty limitation. So one 
should resort to other theoretical approaches such as lin- 
ear response theory.0 However, in fact, electronic trans- 
port in solids is equivalent to the wave transmission of 
electrons through a generalized potential barrier, which 
can be associated with a scattering matrix. In measur- 
ing the conductance of a sample, one always corrects 
the sample to some contacts through perfect leadsa. In 
a two-terminal setup (L=left, R=right), the Landauer- 
Biittiker formulae! states that the current I can be ex- 
pressed as a convolution of the transmission probability 
T and the Fermi distribution function /„ (a = L, R), i.e., 
I=^J T{e)[fLie)~fR{e)]de. The conductance 5 in the 

linear response regime is G = ^ J T{e){—^)de. Such 
a formulation seems more appealing since the transport 
properties are encoded in the corresponding transmission 
probability, which can be calculated by various methods. 

Due to the recent development in nanofabrica- 
tion and material growth technologies, several kinds 
of mesoscopic hybrid structures have been realized 
experimentally. These nanoscale structures include 
normal- metal/superconductor nanostructures (N/S)EI, 
supercoBductor-insulator/superconductor junctions 
(S/I/S)S, superconductor /quantum-dot / superconductor 



transistors (S/QD/S)u, normal metal/superconductisg 
quantum-dot/normal metal transistors (NZSQD/N)ll3, 
ferromagnet/superconductor (F/S) contactEll, supercon- 
ductor/fetromagnet/superconductor sandwich structure 
(S/F/S)E3 and normal-metal/ferromagneti£=(quantum- 
dot/normal metal (N/FQD/N) transistorsEj. In the 
presence |-of a superconductor component, Andreev 
reflectionlU dominates the transport process in the 
low bias case. The imbalance between the spin-up 
and spin-down density of states at the Fermi level for 
ferromagnetic materials introduces the spin-dependent 
transportEj. The conductance of a F/S junction is 
predicted to be smaller or larger than the N/S case, de- 
pending on the ijatio between the spin-up and spin-down 
density of statestS. 

The ngnequilibrium Green's function (NEGF) 
approachtZI has proven to be a powerful technique to 
investigate the transport problem in mesoscopic systems. 
Using the NEGF method, Meir and WingreenEJ have 
derived a Landauer-Biittiker-type formula for transport 
through an interacting normal metal attached to two 
normal leads. Later, some groups have applied ths. 
NEGF method to the N/NQD/SEI and S/NQD/SEI 
cases. In this paper, we extend the NEGF theory 
to a mesoscopic hybrid F/N/S structure, obtaining a 
Landauer-Biittiker-type current formula, which allows us 
to investigate the spin-dependent current and Andreev 
reflecting current in a unified way. 



II. FORMULATION OF THE PROBLEM AND 
CURRENT FORMULAS 

We consider electron tunneling through a ferromag- 
net/normal metal(semiconductor)/superconductor hy- 
brid structure. Under the mean-field approximation, 
the ferromagnet is characterized by a molecular mag- 
netic moHient M, making at an angle 9 with the F/N 
interfaced, while the BCS Hamiltonian is adopted for 
the superconductpt, with an order parameter A describ- 
ing its energy gapE3. In the central region which contains 
a normal metal, we take into consideration various kinds 
of couplings, such as the electron-electron interaction, 
electron-phonon interaction, etc. Then the hamiltonian 
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can be written as 

H = Hp + Hs + -ff_D 

in which 
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(4) 



(5) 



are the hamiltonian for ferromagnet, superconductor, 
central normal metal and tunneling hamiltonian, respec- 
tively. In Eqs. (2-5), /(/^), s{s~^) and c(c+) represent 
the electron annihilation (creation) operator of the ferro- 
magnet, superconductor and normal metal, respectively; 
j^L/R (jgj^otes the tunneling matrix between the normal 
metal and the left ferromagnet (right superconductor). 
Hint is the interaction term in the central normal metal, 
which allows various kinds of interaction to be included. 
Throughout this paper, we assume A is real for conve- 
nience. 

The current flowing into the central region from the left 
ferromagnet can be evaluated from the time evolution of 
the total electron number operator in the left leadll3 



Jl 



dt 



< [Nl,H] > 
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Here we have expressed various kinds of Green's func- 
tions in a generalized Nambu representation, spanned in 
the spin-dependent particle-hole space, in which the spin 
effect and Andreev reflection are considered on the same 
footing and treated in a unified way. The Green's func- 
tions in the generalized Nambu space are of the form 



= ±^E<Y](0/X.(t)(g) 
X,(t)/Y](t') >, 

= T*E^(±^Ti')<[X.(0(8)Y](O 



(7) 



Yt(t')(R)X,(t)>, 



(8) 



in which the four component of the spin-dependent 
particle-hole vector reads 

X,-(X,| Xl X,^ Xl^)\ (9) 

The tunneling matrix takes the form in the generalized 
Nambu space 



kn I p 



(Tkn--\lTpn:] 




V 



^kn-il -^pn-l 














(10) 



kn:] / pn\] ' 



With the help of Dyson's equation, G^^ can be decou- the generalized Nambu representation 
pled into the product of the unperturbed Green's func- 
tion gfefe of the ferromagnet and the Green's function ^ ^ ^ , 
Gran of the Central normal metal regime in the presence 

r^(c.) = 27rT^T^i(c.)TL-*[SI(c^)-S2(^)](12) 
of elastic coupling to the outside world ffl(^) ^ 27rt^,tpfl,('^)T^m = *[Sfl(w) - S?j(c^)](13) 

G^Ji.O = E / rfii[gfc<fc(i,ti)TLG;;.„(ti,t')+gL(t,ii) 

m where (see Appendix) 

Tfc,„G<„(ti,t')]. (11) 
For convenience we introduce the linewidth matrices in Yi^l"' (u) = Tl^^^j ^ (iSf^^^ (14) 
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= ps(w)rfl, 
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with 
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In the above equations, pL| and pLi are the spin- up and 
spin-down density of states in the ferromagnetic region, 
respectively, pg is the normal density of states when the 
superconductor order parameter A = 0, while ps is the 
dimensionless BCS density of states of the superconduc- 
tor. 

Substituting Eq. (11) into (6), we get the following 
compact form for the current 



Jl 



le 



2 — 1,3 



M[G< +fL(GL-G?J]r^}... (25) 



In deriving Eq. (25), we have used the fluctuation- 
dissipation theorem = f (g° — g*^), where f is the Fermi 
distribution function matrix in the generalized Nambu 
representation, given by 



^L/R 



L/R-ij 



(26) 



Here fL/R;ij{x) = Sijlb/Ri!-^ + (-l)VL/fl), with ^l/r 
the chemical potential of the left ferromagnet /right su- 
perconductor. Similarly, one can derive the current flows 
into the right superconductor 
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= f E / ^^{[<^< -f ffl(GL - G?J]f ^}... (27) 



Combining the Eqs. (26) and (27), we get the following 
current formula( J = {Jl — Jr,)/2) 

^ = /M(f^fL"f«fe)(G^,-G^J + 

(f^-f^)G<k, (28) 



i=l,3 • 



which is the central result of this work. It is the analogy 
of Meir-wingreenlla formula in the presence of both ferro- 
magnetic and superconducting leads, providing a frame- 
work to study transport in mesoscopic hybrid interacting 
structures. 

In the absence of interactions within the central region, 
one has the following Keldysh equation 



and the equality 

g: 



a 



i(fLr^ + fflf «), 
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(29) 
(30) 

(31) 
(32) 



With the help of Eqs. (29-32), Eq. (28) can be sim- 
plifled into the following Landauer-Biittiker-type current 
formula for the noninteracting F/N/S mesoscopic hybrid 
structure (setting fij^ — 0, p^ — eV due to gauge invari- 
ance) 



J = J 



NC 



(33) 



in which 
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(35) 



From the expression (33) for the current, one sees that 
the current comes from different physical process: (1) 
jNC includes (a) conventional electron tunneling cur- 
rent, (b) tunneling current of an electron in the ferro- 
magnet into the superconductor either as a hole or as 
a cooper-pair by picking up another electron; (2) J'^ is 
the Andreev reflection current, representing an incident 
electron from the ferromagnet is reflected as a hole back- 
wards into the ferromagnet, with a cooper-pair left in the 
superconductor. At zero temperature, when eV < A, 
J^*-^ becomes zero {ps in is zero when eV < A) 
and the Andreev reflection dominates the current. Ob- 
viously, the current has a strong dependence on the po- 
larization V — '^'"^ , ^"^^ of the ferromagnet and the direc- 
tion 9 of the magnetic moment M, through the coupling 
matrix (notice that the Green's functions of the cen- 
tral region is also implicitly dependent on F^). Equation 
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(28) and Equations (33-35) can be applied to the inter- 
acting and noninteracting mesoscopic F/N/S, N/N/S, 
F/N/N and N/N/N structures, respectively. 

III. EXAMPLES 

In order to illustrate how formulas (33-35) are applied 
to the F/N/S, N/N/S, F/N/N and N/N/N structures, 

I 



From Dyson's equation, G'^/'^ = g''/'^ -1- g^/^'E'^/^G'^/^ , 
one has for Gcc" the expression 

Git = {adj{k)„)det(k)-\ (37) 

where adj{X) denotes the adjoint of X. Note that Gcc" 
is symmetric , i.e., G^^."^- = G^^.^j, as seen from its defi- 
nition or Eq. (36). 

A. current through a F/N/S structure 

For simplicity, we just consider here = 0. Then the 
self-energy matrix S^^" as well as the linewidth matrix 

I 



with 

Bi = detiGl'"), 
B2 = detiGl^"). 

I 



we disciiss in what follows a simple case, in which the 
interaction in the central region is omitted, and only one 
level is relevant to the transport. 

Consider the hamiltonian of the central region, He = 
EcC^c. The unperturbed retarded/advanced Green's func- 
tion gcc" in the 4x4 spin-dependent particle-hole space 
is then 



(36) 



I 

is diagonal, and Gcc" becomes 



where 

(39) 

(40) 

I 

Then the current tunneling through the F/N/S structure 
is 



u! + ec±iO+ 
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i=3A 



dcuihiiu - eV) - hioj + el/)][|GL;i2l' + |GL;34l' 



(41) 
(42) 



The physical imphcations of Equations (41) and (42) 
are apparent. is directly proportional to the spin- 

dependent tunneling matrix due to the coupling be- 
tween the central region and the ferromagnet, and Tj^ as 
well as ps for the superconductor. The Andreev current 
is proportional to the spin-up and spin-down tunneling 
matrix T^^ and F^^, which means that an up(down)- 
spin electron(hole) incident from the ferromagnet will be 
reflected in the interface of N/S and re-enter the ferro- 
magnet as a down(up)-spin hole(electron). In contrast 
to the N/N/S structure (below), the Andreev current 
is strongly dependent on the polarization of the ferro- 
magnet V = P^'^-P^i via the factor rLtrLi . When the 
ferromagnet is fully polarized, i.e., V — I or V = ~l, the 
Andreev current will be zero since no state is available 
for the reflected particle with reverse spin. 



B. current through a N/N/S structure 

When the molecule magnetic moment M is set to 
zero, the ferromagnet becomes a normal metal struc- 
ture. The spin-up and spin-down density of states 
and pli will be the same. Then = T^i — Tl, and 



Gr f^T 
cc:33j "-^c 



Gr 



'^cc;34 — ~^cc;43- Therefore the current through the 



N/N/S system turns into 
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2eri^TLPs 



dL,[fL{u;-eV)-fR{Lo)] 



(lE \G:,af-'2-Re{G:,aiG7c-.i2)],m 



J"" = 



2eT^ f 

J du;[fL{u; - eV) ~fL{io + eV)] 



(44) 



which is the same as the results for the N/NjQD/S hy- 
brid structure derived by Sun and co-workeralj and has 
been studied in detail. 



C. current through a F/N/N structure 

When the order parameter A = 0, the superconductor 
can be viewed as a normal metal. In this case, one can 
observe that ps — ^ and all the non-diagonal elements of 
the full Green's function Gcc'^ of the central region are 
zero. Then the Andreev tunneling current = 0, and 



J is simplified as 



J = J 



NC 



^ J d^[Mco-eV)-MLo)][Tr^^\G:,.^,,\' 
+rLi|GL;33n- (45) 



Since, for normal metal, one can write Tr = (Tri + 
riii)/2, then one gets the following expression for the 
current 
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2h 

eTLjTRi 

2h 

2h 

eT^Li^R-l 



2h 



\G 



cc;33 I 



dujifUuj ~ eV) ~ fR{u;)]\G:,.^,,\^ 

dc.[/L(a>-ey)-/fl(c.)]|GL;33l', 
du;[fL{u;~eV)-fRiu;)]\G:,.^,,\^, 
dLo[fL{to ~ eV) - fRiu)] 

(46) 



The physical meaning of Eq. (46) is very clear. The cur- 
rent into the right normal region comes from two kinds 
of physical process: one is that an electron exits from 
the ferromagnet and enters the central normal region, 
and without spin-flip it tunnels into the right normal 
metal; the other is the spin-flipped tunneling process, 
the up(down)-spin electron leaving the ferromagnet into 
the central normal metal, and after spin-flip, it gets into 
the right normal region with down(up) spin. 



D. current through a N/N/N structure 

Taking advantage of the results for an F/N/N struc- 
ture, if one further assumes F^t — ^li, the F/N/N 
structure then turns into the N/N/N structure, and oue. 
easily recovers from Eq. (46) the well-known expressionc2l 



J = J 



NC 



2eFLFfl, 



dL^ifU^ - eV) ~ fRiLo)]\G:,..,\ 



2e 



diu[h{oJ - eV) - fR{uj)] 

TlTr 

(u;-ec)2 + (Fi+Fj^)2/4' 



(47) 



One sees from Eq. (47) that the transmission probability 



(a;-£,)^+(ri,+r„)V4 



is of the Breit-Wigner form. 
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IV. CONCLUDING REMARKS 

We have given a general Landauer-Biittiker current 
formula, which can be applied to the F/N/S, N/N/S, 
F/N/N and N/N/N structures, even in the presence 
of interactions within the central region. However, one 
finds that one usually can not divide the current into the 
spin-up and spin-down current parts. Then the polar- 
ization in the presence of a ferromagnetic and a super- 
conducting component can not be defined by the current 
polarization. The reason is that the spin-up/down cur- 
rent can not be conserved during the transport process, 
due to the spin flip in the central region and Andreev 
reflection at the N/S interface. For example, a spin-up 
electron from the ferromagnet tunnels through the bar- 
rier into the central region, while another electron with 
spin-down will probably enter the superconductor. It is 
obvious in such case that the spin-up/down current is not 
conserved. 
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Appendix 

In this Appendix we derive the retarded(advanced) self 
energy matrices S^^" and S^^", due to the couphng of 
electrons in the central normal metal to the left ferromag- 
net and the right superconductor. The key step is to get 
first the retarded(advanced) Green's function for the fer- 
romagnet and the superconductor: g^^" and gls"- The 
most convenient way is to diagonalize the hamiltonian 
of the ferromagnet and superconductor. We first calcu- 
late the Green's function for the ferromagnet. Aeplying 
the following Bogoliubov-Valatin transformationcJ to the 
ferromagnet hamiltonian Hp 



fka = cos{e/2)aka - <T sin(0/2)a _ , (48) 



one has 



(49) 



Defining the following retarded(advanced) Green's func- 
tion 

k 

= TiO{±t T t') e-*(''=''+'^*^)(*-*')/f50) 
fe 

r I a 



Zlit't') = T^0i±tTt')Y.<ak4t),alit')> 



(eka+<yM)(t-t')/n 



,(51) 



one finally gets for the ferromagnetic lead the retarded 
(advanced) Green's functions 

- - E ~ = ^anlnnit ~ t')l (52) 



where 

r /a 



g/j.^-^{t-t) = cos ^5ai;e(^"0 



sin^k(re(*-0, (53) 



5//;22(*-0 = cos^^5^i;h(*^0 + 



sin^k&(i-0, (54) 



9ff,zS-^) = cos ^5al;e(^"0 + 



sin^^3:£(i-0, (55) 



g/f.^^{t^t) = cos ^9a\;h{t-t) 



«in^k&(i-0, (56) 



r/a ,, ,/x _ .r/a 







sin 6* r/a 
2 (5'aT;e(^- 


- - 9t 


9tA^-^) 




sin 6* rla u 
2 (5aT:/i(^ 


- - 9:{ 


-r/a ,/N 

g/f.ii{t~t ) 


= .9//;2l(* 


-r/a ,/s 
.9//;23(*-0 


= 0, 


-r/a ,/-, 

5//;34(*-0 


-r/a /, 


-r/a 1, 

9ff-Mi-t) 


= 0. 



(59) 

') 

(60) 

The sum over k can be transformed into an integral, i.e., 
Sfc ~^ I depLa-- Then one arrives at 



.9;;^ii(t - t') = TiSit - t'){cos^ -pL^ + 
9)%2{i-i') = Ti6{t~t'){cos''^-pLi + 

■20 . 

sm ^PL\), 
5?/?33(^-i') = Ti6{t-t'){cos''^-pL^ + 

■20 , 

f/fMit-t') = TiSit-t'){cos^^pL^ + 

■20 ^ 



(61) 



(62) 



(63) 



(64) 
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9ff;lS\^ ~ ^ J — 9ff;SlV' ) — 9ff;24\^ ~ ^ ) 

= 9}%2it-t')=TiS{t-t') 



sin 6 



(PL-f - Pli), 



(65) 



= 9f%S-^')-^^ (66) 
9ffMt -i)= 9ff,Mit - O = 5//;4i(* - * ) 

= 9/;A3it-t')=0- (67) 

The matrix product Lj,,jgjj"(t — t')'Lkrn after 
Fourier transformation J dcue"^'' F{t) yields the re- 

tarded(advanced) self-energy matrix S2^"(a;) Eq. 
(14). 

The procedure to get the retarded(advanced) Green's 
functions for the superconductor lead is similar. After 
the following Bogoliubov-Valatin transformation 



with 



s-pl = P-ph-pi - J^p^l^-, 



l4 = 1^^ + 72^ 



1 



(68) 
(69) 

(70) 
(71) 



the superconductor hamiltonian Hs is diagonalized in the 
following form 



p pa 

2A^ lipi^p. (72) 
p 

Using the notation 

9li-eit - t') = T^O{±t T i') E < bp-it)'blit') > 

p 

p 

9i\i^ - = TiB{±t T t') E < ^p<^(*') > 

p 

k 

one can express the retarded (advanced) Green's function 
matrix for the superconductor lead as 

Elicit - i') = E - = iaU^mnit - t% (75) 



where 



9:i%2{t-t') 

9ss;33it-t) 

^r/a ,i\ 
9ss;44ii-t) 



2 r/a 1 1 
2 r/ 



lit 



'•r/a I, ,i\ 

9ss;12{t-t) 



9ss;34{t - t ) 



'■■r/a /, ,/\ 

9si;13it-t) 



'r/a /, ,/N 
9ss;3lit-t) 



2 r/ 
2 r/ 

Pp9b\ 

'r/a 

9ss-2l\t- 

r/a 
-9bV,eit- 

'r/a /, 

9ss-A3Kt - 

r/a /, 

-9bV,hit - 

^r/a (, 
9ss;14{t - 

— 9ss;24{t — 

'r/a /, 
9ss;32{t- 

9ss;42{t ' 



t') + ^p9lllit-n (76) 

t')+'^yb{:eit-n (77) 

-0 + -XT;l(i-t'), (78) 
-i') + ^p%te(^-*'), (79) 
t') = Ppi^p{9li^{t - t') 
-t')), (80) 
t') = Ppi^pigll^eit - 1') 
-t% (81) 

t')=9:iut-t') 

t') = 0, (82) 

t')=9:i:4i{t-t') 

(83) 



t') 



0. 



Since 



= \ r dep{e-'V^^r/n ^ 

^ J-oo 



I 



Ve2 - A2 

/OO 
de 
-OO 



e\0i\e\-A) 



Ve2 - A2 



(84) 



and 



J dcpfip Up (e 



= - T 

2 7-, 



dCr, 



e2 + A2 



- A2 ^ ^ 

^ Je|g(|6|-A) A^_,,,^, 



Ve^ - A2 e 



(85) 



then one obtains after changing the sum into an in- 
tegral / dcpPg 
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9ss;22it -t)= 9ss;33it " O = 5«s;44(* " * ) 



de 



(86) 
(87) 



Tie{±tTt') J 



6(t-f')/n 



9:i%2it - t') = C{t - t') = gliUt - t') = 0. 



(89) 

(90) 
(91) 



Hence the retarded (advanced) self-energy matrix transformation J da;e*"*F(t). 
S^''"(a;) Eq. (15) is similarly obtained from the di- 
rect matrix product R|„gls°(t — t')'Rpm after Fourier 
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